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In this letter we show that although the application of standard Lyapunov analysis predicts that completely 
integrable Kepler motion is unstable, the geometrical analysis of Horwitz et al [5] predicts the observed 
stability. This seems to us to provide evidence for both the incompleteness of the standard Lyapunov analysis 
and the strength of the geometrical analysis. Moreover, we apply this approach to the three body problem in 
which the third body is restricted to move on a circle of large radius which induces an adiabatic time dependent 
potential on the second body. This causes the second body to move in a very interesting and intricate but 
periodic trajectory; however, the standard Lyapunov analysis, as well as methods based on the parametric 
variation of curvature associated with the Jacobi metric, incorrectly predict chaotic behavior. The geometric 
approach predicts the correct stable motion in this case as well. 



PACS numbers: 45.20.Jj, 47.10.Df, 05.45.-a, 05.45.Gg 



Several relatively recent papers have used a geometric 
approach to describe Hamiltonian chaos, implementing tools 
from Riemannian geometry [1-3]. The natural motions of 
Hamiltonian systems are viewed as geodesies on a 
Riemannian space associated with a metric g (often 
associated with cither the Jacobi or the Eisenhart metric). 
Stable motions are thus de.ned by the ctirvature properties of 
the manifold. Casetti et al [2] derived from the Jacobi.Levi- 
Civita (JLC) equation an effective stability equation which 
formally describes a stochastic oscillator. They conjectured 
that some 'average' global geometric property should provide 
information about the degree of chaos and applied the 
geometric method to the Fermi-Pasta-Ulam beta-model and 
to a chain of coupled oscillators. Safaai and Saadat [4] have 
used the geometric method developed in [1-3] in an attempt 
to predict chaos in the restricted three body problem. They 
conclude that the .fluctuations of the curvature of the 
manifold along the geodesies yield parametric instability of 
the trajectories and cause chaos. They also calculated a 
positive Lyapunov exponent. It has recently been shown [5] 
that there is a possibility to characterize instability in 
Hamiltonian systems by a geometrical approach which takes 
its point of origin in the curvature associated with a 
conformal Riemann metric tensor which is apphcable to a 
large class of potential models. This approach appears to be 
more sensitive than computing Lyapunov exponents or the 
use of the Jacobi metric. This enables one to associate 
instability with a negative dynamical curvature appropriate 
for the geodesic motion different from that defined by the 
Jacobi mefric, for which the line element is proportional to 
the time x and not, as for the Jacobi metric, the action. This 
establishes a natural connection between chaotic Hainiltonian 
flows and Anosov flows [6]. 

The geomettical method that we shall use [5] studies the 
stability of a Hamiltonian system of the form (we use the 
summation convention) 

H = ^ + V(x) (1) 
2M ^ ' 

where V is a ftmction of space variables alone, by inttoducing 
a second Hamiltonian of the type considered by Gutzwiller 
[7, 8] 



H = —gup'p' 
2M 



(2) 



where gjj is a function of the coordinates alone. 
Hamilton's equations applied to Eq. (2) result in 

X = —Ti Xj^X^ ; (3) 

where 



2 dx„ dx^ oxt 



-} 



(4) 



and g'^ is the inverse of gy . 
Horwitz, et al [5] take g^yto be of conformal form 



Sij 



E 



(5) 



E-V{x) ^' 

on the hypersurface defined by H = E = const, 
resuhing in a formal equivalence between the two 
Hamiltonians on the same energy shell E . 
The velocity field defined by 



U-' =X^ = g'^Xj = 



M 



satisfies the geodesic equation 
mn 

where 

1 Ik SSmn 



mn r> S 



(6) 



(7) 



(8) 



It is shown in [5] that in the coordinates for which 
gjj has the form (5), (8) coincides with the Hamilton 

equations obtained irom the Hamiltonian form (1), and 
thus provides a geometrical embedding of this 



Hamiltonian motion. The geodesic deviation E,^ 
between two such orbits x'(?)and satisfies 



M J,/ 



DMt' 



(9) 



where the covariant derivative of ^' is defined by 



(10) 



M 



is the corresponding covariant derivative along 



the geodesic curve, and 

dML dMl 
R = - 



qm 



qn 



Mqmn 



dx" dx" 



.(11) 

With the conforiml metric in noncovariant form (5), 
the geodesic deviation equation (9) becomes 



M 



where the matrix V is given by 

V 3 dV dV ^ \ d^V 
Mv^ dx' dx' M dx'dx' 



and 



(12) 



(13) 



(14) 



with =x\ defining a projection into a direction 

orthogonal to v' . 

It was argued in [5] that instability should occur if at 
least one of the eigenvalues of V is negative. This 
criterion was found to be effective for many examples 
[5]. 

The analysis given above was carried out with a time 
independent potential function V . We now study the 
consequences of allowing a weak time dependence in 
V . In this case E is not precisely conserved, but we 
shall make an adiabatic approximation in which we 
consider Eio be time independent. The time 
dependence of V appears in the form of the geodesies 
generated by the Gutzwiller Hamiltonian (2) as an 
additional term with partial time derivative of the 
metric, but in the geometrical imbedding of the 
Hamiltonian motion, the geodesic evolution of the 
velocity vector field has the same form (7) as for the 
time independent potential problem. The formulation 
of geodesic deviation in terms of the second covariant 
derivative, however, introduces another time 
derivative, introducing an additional term in the 
stability matrix. Following the methods applied above, 
one finds that the geodesic deviation now obeys the 
relation 



(15) 



which is just the result obtained in [5], with the 
additional term 



dM 



rt -t i^r _^ ^ Srt -t nr 



dt 



2 dtdx. 



(16) 



Using the conformal form of gjj with time dependent 
potential, the right hand side contains 



2E 



dV dV 



+ - 



E 



dtdXg {E-V{x,t)f dt dx^ ' {E-V{x,t)f dx^dt 

(17) 

Near the boundary of the physical region, these terms 
can become large, but they are generally small away 
from these boundaries since they involve the time 
derivative of F , in the same approximation in which 
E is constant (the consistency of this approximation 
will be discussed elsewhere). In our computations here 
these conditions hold, and we therefore compute 
stability using the same formulas as for the time 
independent problem. 

We now consider the following Lagrangian of a 
restricted three body system: 



L = T-V, T = -m 



dt 



1 2 2 

V-—m=rf 0); 
2 



2 



2 

4;r m, 



47r m^nij 



'12 



.(18) 

In this system three bodies are involved: a large mass 
M (the "Sun") assumed to be set at the origin (see 
Fig.l), a mass ("Earth") which its frajectory we 

wish to study, located at r , and a third mass ntj (a 
"Jupiter") located at Pj circling the sun with angular 
frequency CDj such that 

rj = rj [cos[cOjt),sm[cOjt)) , (19) 
Vj is assumed to be constant. 




FIG. 1: Depiction of the restricted three body problem 





FIG. 2: The Kepler orbit [ntj = o) 

Finally is the distance between the "Earth" and 
"Jupiter". The following parameters are assumed: 



rj= 5.203, 0)j 



nijQ = 317.9 -/Mg 



In: 



=2.96-10" 



Two cases are studied: in the first case we ignore 
Jupiter, Wy = ; (the classical Kepler problem) and in 

a second case we assume a "large Jupiter" 
rrij — lOOOm^o . The trajectories of the earth in those 
two cases are depicted in Figs. 2 and 3 for the initial 
conditions F(o)=(1,0), /"^(o) = (0,1) . 




FIG. 3: The earth orbit for the restricted three body problem 
with nij = 1000?Wyo during 100 time units 



In the geometrical approach we write the matrix V ; 
Eq.(13), as: 



kl- 



1 

m„ 



dV dV 



2[E -V) dXj^ dxi dxi^dxj 



(21) 

where E is now assumed to be adiabatically constant. 

For the Kepler case, one eigenvalue is always 
positive; the second eigenvalue is depicted in Fig. 4. 
Clearly it becomes negative during certain short time 
intervals. We found that this occurs at the perihelia of 
the orbits. A local Lyapunov instability is not a cause 
for concern since the orbit may be globally stable. 
However, Fig. 4 reveals that the Lyapunov eigenvalue 
is indeed negative along the entire orbit although its 
negative value is less pronounced between perihelia 
than at the perihelia. Thus according to Lyapunov 
analysis the Kepler problem is unstable and could 
become chaotic. This is clearly a false result since the 
Kepler problem is integrable. 

Yet the "geometric eigenvalue", depicted in Fig.5, is 
found to be positive in between perihelia, hence the 
trajectory is predicted to be globally stable. The 
trajectory only spends a short time at the unstable 
perihelia and thus instability does not have enough 
time to take place. This is a clear indication that the 
geometrical analysis is preferable. 

The trajectory of the restricted three body problem is 
shown in Fig.3. This trajectory is quite complicated, 
however, it appears to be periodic rather than chaotic. 
The Lyapunov eigenvalues and the geometric 
eigenvalues are depicted in Fig.6 and Fig. 7. 

The eigenvalues in the restricted three body problem 
have a different behaviour than those of the Kepler 
problem. Nevertheless, the same conclusions are 
obtained: Lyapunov analysis predicts instability and 
possible chaos, geometrical analysis predicts stability. 
The application by Safaai et al [4] of the methods of 
Cassetti and Pettini also appear to predict chaotic 
behaviour, which is not in agreement with our 
simulations. 



We now study the stability of each of the trajectories 
described in the previous section using the standard 
Lyapunov analysis and evaluating the eigenvalues of V 
in Eq.(13). The eigenvalues of this matrix are supposed 
to determine the stability of the trajectory. A negative 
eigenvalue indicates an unstable trajectory while a 
positive eigenvalue indicates stability. It can be easily 
be seen that the Lyapunov exponent depends on the 
eigenvalues of the stability matrix (denoted by 



'Lyapunov 



V, 



in the figures) 



1 



X , X2 = y, 



(20) 

the square of the Lyapunov eigenvalues are given by 
the negative of the eigenvalues of this matrix. 
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FIG. 4: The eigenvalues of the Lyapunov stabihty matrix 
(Eq.(20)) for Kepler problem, nij = . Upper left: during 

Itime unit; upper right: passing a perihelion; middle: between 
perihelia. 



FIG. 6: The eigenvalues of the Lyapunov stability matrix 
(Eq.(20)) for the restricted three body problem with 
my =1000m|o- Upper left: during 1 time unit; upper right: 

passing a perihelion; middle: between perihelia. 
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FIG. 5: The eigenvalues of the "Geometry" stability matrix 
(Eq.(21)) for Kepler problem, nij = . Upper left: during 1 

time unit; upper right: passing a perihelion; middle: between 
perihelia. 



FIG. 7: The eigenvalues of the "Geometry" stability matrix 
(Eq.(21)) for the restricted three body problem with 
nij = lOOOmyo . Upper left: during 1 time unit; upper right: 

pEissing a perihelion; middle: between perihelia. 
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